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We prove that Nambu-Jona-Lasinio model is an exact description of infrared QCD deriving it 
from QCD Lagrangian. The model we obtain is renormalizable and confining but, taking very small 
momenta fixes completely all the parameters of the Nambu-Jona-Lasinio model through those of 
QCD. The choice of the infrared propagator is done consistently with recent numerical results from 
lattice and Dyson- Schwinger equations for Yang-Mills theory. The model we get coincides, once 
the ultraviolet contribution is removed, with the one proposed by Langfeld, Kettner and Reinhardt 
[Nucl. Phys. A 608, 331 (1996)]. 
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©J). I. INTRODUCTION 

' Studies of low energy strong interaction processes generally meet severe difficulties due to the known impossibility to 
, manage the right theory, i.e. Quantum ChromoDynamics (QCD), in this limit. So, a lot of phenomenological models 
fNj ■ have been devised to cope with this problem. One of the most used is the Nambu-Jona-Lasinio model 0, 0, H, BL IE El • 
Despite some severe difficulties this model has, a large body of literature developed about as it reproduces most of 
the observed patterns in nuclear physics. The main ones are that this model is not renormalizable, one has contact 
interaction, and fails to maintain quarks confined. 

Some years ago, Reinhardt, Kettner and Langfeld provided a renormalizable Nambu-Jona-Lasinio model once 
the one-gluon exchange propagator is known. At the time this latter information was lacking but the authors put 
, out a plausible model for this propagator and were able to perform significant computations in nuclear physics. This 
model, as we will show, implies that a strong similarity exists, at least in the infrared, between QCD and Quantum 
ElectroDynamics (QED). Indeed, this deep insight is correct as we will see implying a justification for authors' 
computations and also for a lot of published literature using Nambu-Jona-Lasinio model that represents the right 
^\ ' behavior of hadronic matter in the low energy limit. 

In order to have a consistent argument, two main points should be addressed. Firstly, one should know if a Green 
function method does hold also in a strong coupling regime. We have proved this in a series of papers [1, 0, 
C "— ; , showing that exists an expansion for nonlinear equations, that holds at small times, implying the Green function of 
these equations. When a Fourier transform is done, one gets powers of the energy giving a series that holds when 
energy goes to zero as we need. Secondly, one has to have a proper gluon propagator, that is the right one from 
Yang-Mills theory. 

About the question of the behavior of the gluon propagator in the infrared there is currently a lot of activity. A 
first understanding about this matter was put forward by Alkofer and von Smekal (a review is that, introducing 
a truncation to the Dyson-Schwinger equations for a pure Yang-Mills, obtained a gluon propagator going to zero in 
the infrared. A older work due to Cornwall [l2| has showed a contradictory result with a gluon propagator reaching 
a constant value as the energy goes to zero. Both scenarios implied a running coupling reaching a fixed point in the 
infrared. The present situation has reached a different grasp of the situation mostly through lattice computations 
and numerical solution of Dyson-Schwinger equations. Recent lattice computations proved that indeed the gluon 
propagator reaches a non-zero constant value in the infrared but the running coupling does not reach any fixed value 
but rather it goes to zero [HI, 0, [Hj]- Same happens in the numerical solution of the Dyson-Schwinger equations 
G3] • A running coupli ng g oing to zero has been also seen on lattice computations in [To ] and with an analysis of 
experimental data p^ . [20l . |2l| . So, consensus is going to be reached about a gluon propagator going to a non-zero 
value as momentum goes to zero, the ghost propagator being the same of that of a free particle, and a running 
coupling, also with different definitions, going to zero. This scenario was already well described by Boucaud's group 
with lattice computatio ns [22l [23l |24| . This emerging scenario is in complete agreement with the one we have put out 
recently [25L l26l. WA l28l [2~°| based on Q. This gives a gluon propagator in a closed form in very good agreement, as 
we will show, with lattice computations and with the numerical solution of the Dyson-Schwinger equations while the 
running coupling obtained with Callan-Symanzik equation [29| is in agreement with lattice computations [l8l ]. 
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The paper is so structured. In [IT] we outline our argument with a reformulation of QED. In IIIII we apply the same 
argument to QCD obtaining a rcnormalizable Nambu-Jona-Lasinio model to be applied in the infrared. In IIVI we 
discuss the question of the ghost field and Gribov copies in the infrared limit where both seem to be excluded. In 
IVl we show how our choice of the gluon propagator is indeed a consistent one with respect to lattice and numerical 
results. In I VII we obtain the contact interaction proper to the original formulation of the Nambu-Jona-Lasinio model 
fixing all the parameters of the theory. We give also some discussion of the ground state of QCD with respect to a 
pure Yang-Mills theory. Finally, in I VIII conclusions are given. 



II. QUANTUM ELECTRODYNAMICS 



Let us consider a classical action for the interaction between a charged particle and photons 



S = 



1 



1 



V>(*7 ■ d - ej ■ A - m)tp - 7 (d^A u - d fJ ,A u )(d^A v - d^A") - — (d ■ A) 



2a 



(1) 



This action gives the following equations 



(ij ■ d — e7 • A — m)ip = 



d^d 11 A" 



1 - - ) ff'&'Au 



The solution for the electromagnetic field can be obtained immediately through Green function as 

A^x) = ej dS^ix - y)4>(yhMv) 

being 



if - (1 - a) 



k^k v 
k 2 
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the photon propagator. Substitution of eq.([3]) into action ((T|) finally gives 

1 



S 



ip(x)(ij ■ d - m)ip(x) - -e ip(x)^^(x) I d yD^(x - y)ip(y)juip(y) 



(2) 



(3) 



(4) 



(5) 



where we can see the appearance of a quartic term. At this stage we can do quantum field theory with both actions 
and so we are granted that both are renormalizable. But we can also see that for QED, due to the photon being 
massless, there is not an immediate way to get a contact quartic interaction out of this latter action. Things can 
be quite different in QCD as we know that this theory in the infrared develops a mass gap. Finally we note that a 
generic perturbation term with a propagator in agreement with the Kallcn-Lchman form 



p 2 — + ie 



(0) 



does not spoil renormalizability of the theory being a sum of renormalizable terms. Indeed, it is enough to have an 
infrared gluon propagator reaching a non-null value at zero momenta to recover a Nambu-Jona-Lasinio model (e.g. 
see [30] ) - Our aim in the following is to derive such a formulation for QCD that holds in the infrared. As we will 
show, we have all one needs to reach this aim. 



III. QUANTUM CHROMODYNAMICS 

At this stage we work out a classical model. The needs for quantization impose the introduction of a Fadeev-Popov 
ghost and to fix the question of Gribov copies. Indeed, in the infrared limit that interests us here the ghost field is 
proved to decouple from the gluon field and the Gribov copies are harmless. We will discuss all this matter in sec. 
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So, the action for QCD can be written as 
S = Id 



L q ^ 



being q — u,d, . . . the flavor index and 

G^ = d„Ai-a v A a M +gf abB A b tl Al 
the field tensor. This action produces the equations 



(7) 



(8) 



ij-d- g—~i ■ A a -m q )q = 



?r hc f cde A h »A%Al = gY,q^luq 



(9) 



At this stage we use the small time approximation devised in [f| [T(| and write down the solution for the field 
equations as 



A%(x) « g J dSDf v {x - y) J2q(y)¥- 7 « q (y) 
while for the gluon propagator we take the form 



ah 



V»» - (l-a)-p- 



G(k). 



(10) 



(11) 



This approximation implies that we get a theory that holds at small momenta. Indeed, this is just the leading order 
of a small time series [1(|. Then we can write the QCD action that holds in the infrared limit as 



^2q(x){i-f ■ d~ m q )q{x) - -g 2 ^q{x) —l^qix) j d 4 yD a Jl(x - y)q{y)—-f u q'(y) 



(12) 



We just note that Nambu-Jona-Lasinio model is recovered if the gluon propagator has a mass gap m. In this case 
one should have e.g. 



G(k) = 



k 2 — to 2 + ie 



(13) 



being m the mass gap. Then, in the infrared m 3> k and then G{k) constant giving back a propagator proportional 
to 8 4 (x — y) producing the required contact interaction while the mass gap becomes a natural cut-off for the theory. 
As already said, for QCD this is indeed the situation. 

Then, our next aim is to fix the form of G(k) in agreement with the current scenario in lattice computations for a 
pure Yang-Mills theory. 



IV. 



GHOST FIELD AND GRIBOV COPIES 



In order to have an idea on how to fix the gauge in the infrared limit, we need to know some results that have 
been obtained recently on the lattice. It should be said that, in some cases, analytical proofs exist but we will avoid 
these as theoretical results need some time before their correctness is properly assessed. So, we heavily rely on lattice 
computations that presently, with the use of huge volumes as we will see in sec[V] have definitely set the scenario 
about the propagators at lower momenta. 

The question of the importance of Gribov copies in this limit has been present since the start of studies of QCD 
infrared limit. Gribov copies are not important in the ultraviolet limit as otherwise we have had serious problems 
to unveil asymptotic freedom and compute higher order corrections, a well-acquired technology these days [31(. In 
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the infrared limit there has been some theoretical work to understand confinement and this set the agenda for people 
working on the lattice as people have had to understand the relevance of Gribov scenario in this framework. 

The most important works on this matter have been published by Ilgenfritz, Miiller-Preussker and collaborators 
[32l [H| quite recently where they consider different volumes and (3 to understand the effect of Gribov copies on gluon 
and ghost propagators. Their conclusions are somewhat unexpected: They do not see any effect of Gribov copies 
on the gluon propagator but the effect appears quite significant for the ghost propagator. From their later paper 
[33| they show as the effect of the Gribov copies is even smaller as the volume enlarges. Indeed, eigenvalues of the 
Fadeev-Popov operator are seen accumulating to zero while the ghost propagator saturates at smaller volumes where 
the effect of Gribov copies seems more effective. So, the question to be understood after these analyses is how relevant 
is the ghost propagator in the behavior of Yang-Mills theory in the infrared. 

The answer appeared quite recently in a pair of beautiful papers due to Cucchieri and Mendes (341 . |35| . They firstly 
showed that the gluon propagator must reach a finite non-null value in the infrared limit [34| . Then, they proved that 
the ghost propagator behaves as that of a free particle [H[ , that means that it decouples from the gluon field and we 
can avoid to care about it in the infrared limit. This means that Gribov copies are not relevant at lower momenta in 
the same way as happens in the ultraviolet limit. So, we conclude that our starting Lagrangian is well supported by 
lattice computations and no further constraints are needed to fix it. 



V. GLUON PROPAGATOR 

The choice of the gluon propagator must be consistent with recent lattice computations [HI, [lj, [l5[ and in perfect 
agreement with the numerical solution of the Dyson-Schwinger equations (l6| . For our aims we do not consider the 
corresponding numerical solution on a torus [17j |. We just note that in this case a strongly infrared depressed running 
coupling is observed. The next most important point is that the gluon propagator must satisfy Callan-Symanzik 
equations and the corresponding beta function should give a running coupling going to zero in the infrared as is seen 
on lattice. Finally, we do not have to consider coupling with the ghost as this is seen to behave as a free particle and 
so decouples from the gluon field. 

We have proved in a series of a paper [2^, [26|, [U [H, [29| that a pure Yang-Mills theory admits a gluon propagator 
given by 

<W-Eld*TE (14) 

ri=0 ™ 



being 



7T 2 f~l) n + 1 e -("+3) 7r 

g " = ( 2 " +l Wi;e-(»*i>- ' (15) 



m n= [n+D-^r ) n (Hi) 



and 

m_ = I 7? -A- . 

2 J K(i) V 2 

being N the number of colors, fi an integration constant to be fixed experimentally and K(i) the constant 

K(i) = { 2 - d ° « 1.3111028777. (17) 
Jo VI + sin 2 

This propagator is the same of a massless scalar field theory obtained in the same limit in 0. They coincide as the 
only useful solution of a Yang-Mills propagator to build a quantum field theory is the one having all equal components 
[2^ | . The reason for this mapping between a massless scalar field and Yang-Mills theory in the infrared limit relies on 
the fact proved in the eighties that most solutions of Yang- Mills quantum mechanics are generally chaotic [H, [37], [38| 
and so, completely useless to build a quantum field theory out them. So, these prove to be the only classical solutions 
of Yang-Mills theory to be used in the infrared limit. Indeed, these solutions are exact solutions of classical Yang-Mills 
equations producing a mass gap out of massless equations of motion [3t| . In this case the Yang-Mills theory coincide 
with the massless scalar field theory provided one substitutes the coupling A with the 't Hooft coupling Ng 2 . 

This propagator fulfills all the requirements seen on lattice computations. It reaches a finite non-null value at lower 
momenta and recovers the ultraviolet behavior being J^^Lo ^ n ~ — 1- ^ should a l so De noted that 

G(o) = -V- r- — ' - 3.761402959--. (18) 

w cr^2n+l l + e -( 2n + 1 ) ,r a v ; 
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FIG. 1: Gluon propagator (circles) compared to lattice data (stars): a) Q (13.2/m) 4 ; b) (19.2/m) 4 ; c) PJ] (27/m) 4 ; d) 
[T^ | Dyson-Schwinger. 



where we have introduced the string tension as y/a = (TV g 2 \x 4 /2) 4 . Eq. (fT5|) introduces a fundamental constant that 
defines the Nambu-Jona-Lasinio model from QCD. 

In order to give a complete proof of the fact that this is indeed a proper choice for the gluon propagator, we want 
to compare it with lattice results and the numerical solution of the Dyson-Schwinger equation. While in the former 
case we expect some error due to volume effects in the lattice computations that should appear to decrease with 
increasing volume, for the numerical solution of the Dyson-Scwhinger equation we must have a complete coincidence. 
We will show that is indeed what happens. We emphasize that in order to do this comparison we only need to fix the 
parameter /i that gives the gluon mass. But we are comparing non homogeneous data belonging to different research 
groups so the value at zero momenta can be different. This should be expected but the only relevant thing is that 
fixing just the value of our parameter we are able to get excellent agreement with lattice and numerical data. 

For lattice computations we have that in [l| the gluon propagator is computed with the largest volume (27/to) 4 , 
the next one is [HI (I9.2/m) 4 and finally [l4| (I3.2/m) 4 . For the numerical solution of the Dyson-Schwinger equations 
(l6| the curves practically coincide. But these computations have not all the same a being different groups. What 
we want to see is that, independently on the value of a, in any case the fits improve increasing the volume. This is 
exactly what happens as can be seen in fig. [T] where comparisons with all propagators is done. 

Indeed, one can see from fig. [2] that the error is significant in the intermediate energy region and not for infrared 
but this error significantl y d ecreases going from the smallest to the largest volume reaching a fairl y g ood agreement 
with the results given in jl 31 ] and a perfect agreement with numerical Dyson-Schwinger equations 16]. The trend is 
that the lattice curve indeed converges toward the propagator we have given here at increasing volume. 

Then, we can assume that the propagator (JT?)) represents correctly the behavior of Yang-Mills theory in the infrared 
and we can take this into our rewriting of QCD given in the action (| 1 2|) . As this propagator is just a sum of massive 
free propagators, we are able to get a Nambu-Jona-Lasinio model describing low energy QCD. 



VI. NAMBU-JONA-LASINIO MODEL FROM QCD 



At this stage we can write down the action of a renormalizablc Nambu-Jona-Lasinio model as 



Snjl ~ / d 4 x 



r ^ r ya 

J d ^ x ^Z^)^! ■ d ^ m q) ( l{ x ) ~ ^^^—^qix) J d 4 yG(x - y)q' {y)—l^q' {y) 



(19) 
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FIG. 2: Percentage error for our propagator with respect to lattice and numerical data. It is seen how increasing volume makes 
converging the data toward our solution. 



where use has been made of current conservation that grants independence on the choice of the gauge in the propagator. 
Here G(x — y) is just the Fourier transform of the propagator in eq. (fT4]) . We emphasize that this model should holds 
just at lower momenta. This situation is quite different from the case of QED. 

At this stage we just notice that our propagator (|14| implies the exchange of particles with masses given in eq. (fT6|) 
and this are not gluons but glueballs that are not confined but would be stable particles if the ground state of QCD 
would be higher than the lowest of the masses of these particles. This is not the case in the real world and glueballs 
decay. The lowest glueball state has an energy higher than the tt meson mass. Then, 7r mesons should appear as 
decay products of glueballs. With the Nambu-Jona-Lasinio model we should be able to compute such decay rates. 

We point out that model given in eq. (fT!?)) is renormalizable due to its strict analogy with QED (0) , the only difference 
being a sum of massive propagators that in the limit of high momenta reduces to a sum of QED Lagrangians. But we 
can also conjecture that this model is confining due to the nature of Yukawa propagators involved into the interaction 
terms. These properties are not maintained with a contact interaction. Indeed, this can be recovered by taking the 
limit of small momenta (infrared limit k — * 0) giving after a Fourier transform 



1 



G(x-y) w3.76- -6 4 {x-y) 



(20) 



that gives back the original Nambu-Jona-Lasinio model 



S: 



d 4 x 



q(x){iH ■ d - m q )q(x) - -G NJL 2J q(x) —^q{x)^{x) —^^{x) 
9 ~ q,q' 



(21) 



being 



Gnjl ~ 3.76 • — 
a 



(22) 



the Nambu-Jona-Lasinio constant linked to the experimental value of the string tension a. We can see a strong 
similarity with the electro-weak part of the standard model that in the low energy limit recovers the Fermi model. 
Here we get again a non-renormalizable model and we do not have to worry about the missing of confinement as 
at increasing momentum we have to use the more general model given in eq. (|19[) and correct it when we are at the 
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boundary of the intermediate energy region. Besides, the model has a natural cut-off originating from the mass 
spectrum of the Yang-Mills theory. Indeed we have assumed 2 K(i) ^ & an d we can take for the Nambu- Jona- 
Lasinio cut-off 

A NJL = (23) 

that is about the mass of the a particle (or f0(600)) (about 527 MeV) when we take for yfa the value 440 MeV 
currently chosen on lattice computations. 

We emphasize that when |fc| > Ajvjl, our approximation fails and higher order corrections are needed. 

As a final consideration we note that this is exactly the model presented in Ref.[7J] when the ultraviolet part is 
removed. Then, we can see that these authors hit the right model and all their conclusions apply provided we consider 
a momentum |fc| < Ajvjx- 



VII. CONCLUSIONS 



We have shown how, from QCD, a renormalizable Nambu-Jona-Lasinio can be derived. This is possible as, in the 
infrared limit, we can still adopt the concept of a Green function. The existence of a mass gap that changes the 
propagator into a massive one, or a sum of free massive propagators, grants the emerging of a contact interaction in 
the small momenta limit. 

It is really surprising the insight of Nambu and Jona-Lasinio that about fifty years ago were able to postulate a 
model that turns out to be the right one to describe strong interactions in the lower energy limit. A large body of 
literature has been published about and the main conclusion is that this research activity has been well deserved. 
On the same ground one should put the more recent work by Langfeld, Kettner and Rcinhardt that postulated a 
renormalizable version of the Nambu-Jona-Lasinio model turning out exactly the one we got from QCD. 

With this formulation of low energy QCD, the possibility is given to perform all kind of computations one needs 
to understand nuclear matter. The model we obtained relies on some experimental values but these are exactly the 
same needed for QCD as it should be. 
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